library(knitr)
library(rmdformats)
library(ggplot2)
library(ggpubr)
library(GGally)
## Registered S3 method overwritten by 'GGally':
## method from
## +.gg ggplot2
library(car)
## Loading required package: carData
library(tidyverse)
library(lme4)
library(lmerTest)
library("MuMIn")
library(lmtest)
library(boot)
Read dataset
AllSubs_NeuralActivation <- read.csv('/Users/luisalvarez/Documents/GitHub/RM_Thesis_Neuroforecasting/ProcessedData/AllSubs_NeuralActivation_Aggregate_Combined_Comedy_clean.csv')
# Define aggregate variables.
All_Gross_W1_log <- log(AllSubs_NeuralActivation$Gross_US_W1_num)
All_Theaters_W1 <- AllSubs_NeuralActivation$Theaters_US_W1_num
All_Gross_M1_log <- log(AllSubs_NeuralActivation$Gross_US_M1)
All_Theaters_M1 <- AllSubs_NeuralActivation$Theaters_US_M1
# Define affect variables.
All_PA <- AllSubs_NeuralActivation$Pos_arousal_scaled
All_NA <- AllSubs_NeuralActivation$Neg_arousal_scaled
FW_US_M1_df <- data.frame(All_Gross_W1_log, All_Theaters_W1)
FM_US_M14_df <- data.frame(All_Gross_M1_log, All_Theaters_M1)
# Define ISC variables.
All_NAcc_ISC <- AllSubs_NeuralActivation$NAcc_ISC
All_AIns_ISC <- AllSubs_NeuralActivation$AIns_ISC
All_MPFC_ISC <- AllSubs_NeuralActivation$MPFC_ISC
# Define models.
FW_US_M2_df <- data.frame(All_NAcc_ISC, All_AIns_ISC, All_MPFC_ISC)
FW_US_M3_df <- data.frame(All_Gross_W1_log, All_PA, All_NA, All_NAcc_ISC, All_AIns_ISC, All_MPFC_ISC)
FM_US_M16_df <- data.frame(All_Gross_M1_log, All_PA, All_NA, All_NAcc_ISC, All_AIns_ISC, All_MPFC_ISC)
# Define whole variables.
All_NAcc_whole <- AllSubs_NeuralActivation$NAcc_whole
All_AIns_whole <- AllSubs_NeuralActivation$AIns_whole
All_MPFC_whole <- AllSubs_NeuralActivation$MPFC_whole
# Define models.
FW_US_M4_df <- data.frame(All_NAcc_whole, All_AIns_whole, All_MPFC_whole)
FW_US_M5_df <- data.frame(All_Gross_W1_log, All_PA, All_NA, All_NAcc_whole, All_AIns_whole, All_MPFC_whole)
FM_US_M18_df <- data.frame(All_Gross_M1_log, All_PA, All_NA, All_NAcc_whole, All_AIns_whole, All_MPFC_whole)
# Define onset variables.
All_NAcc_onset <- AllSubs_NeuralActivation$NAcc_onset
All_AIns_onset <- AllSubs_NeuralActivation$AIns_onset
All_MPFC_onset <- AllSubs_NeuralActivation$MPFC_onset
# Define models.
FW_US_M6_df <- data.frame(All_NAcc_onset, All_AIns_onset, All_MPFC_onset)
FW_US_M7_df <- data.frame(All_Gross_W1_log, All_PA, All_NA, All_NAcc_onset, All_AIns_onset, All_MPFC_onset)
FM_US_M20_df <- data.frame(All_Gross_M1_log, All_PA, All_NA, All_NAcc_onset, All_AIns_onset, All_MPFC_onset)
# Define middle variables.
All_NAcc_middle <- AllSubs_NeuralActivation$NAcc_middle
All_AIns_middle <- AllSubs_NeuralActivation$AIns_middle
All_MPFC_middle <- AllSubs_NeuralActivation$MPFC_middle
# Define models.
FW_US_M8_df <- data.frame(All_NAcc_middle, All_AIns_middle, All_MPFC_middle)
FW_US_M9_df <- data.frame(All_Gross_W1_log, All_PA, All_NA, All_NAcc_middle, All_AIns_middle, All_MPFC_middle)
FM_US_M22_df <- data.frame(All_Gross_M1_log, All_PA, All_NA, All_NAcc_middle, All_AIns_middle, All_MPFC_middle)
# Define offset variables.
All_NAcc_offset <- AllSubs_NeuralActivation$NAcc_offset
All_AIns_offset <- AllSubs_NeuralActivation$AIns_offset
All_MPFC_offset <- AllSubs_NeuralActivation$MPFC_offset
# Define models.
FW_US_M10_df <- data.frame(All_NAcc_offset, All_AIns_offset, All_MPFC_offset)
FW_US_M11_df <- data.frame(All_Gross_W1_log, All_PA, All_NA, All_NAcc_offset, All_AIns_offset, All_MPFC_offset)
FM_US_M24_df <- data.frame(All_Gross_M1_log, All_PA, All_NA, All_NAcc_offset, All_AIns_offset, All_MPFC_offset)
# Seq models.
FW_US_M12_df <- data.frame(All_NAcc_onset, All_AIns_middle, All_MPFC_offset)
FW_US_M13_df <- data.frame(All_Gross_W1_log, All_PA, All_NA, All_NAcc_onset, All_AIns_middle, All_MPFC_offset)
FM_US_M26_df <- data.frame(All_Gross_M1_log, All_PA, All_NA, All_NAcc_onset, All_AIns_middle, All_MPFC_offset)
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ Theaters_US_W1_num + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.21830 -0.37441 0.02446 0.39529 1.19424
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.4675794 0.4474662 27.863 2.83e-12 ***
## Theaters_US_W1_num 0.0011595 0.0001464 7.919 4.17e-06 ***
## scale(Pos_arousal_scaled) -0.0631872 0.1922402 -0.329 0.748
## scale(Neg_arousal_scaled) 0.2691857 0.1889371 1.425 0.180
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7255 on 12 degrees of freedom
## Multiple R-squared: 0.8442, Adjusted R-squared: 0.8052
## F-statistic: 21.67 on 3 and 12 DF, p-value: 3.905e-05
## Warning: 'r.squaredGLMM' now calculates a revised statistic. See the help page.
## R2m R2c
## [1,] 0.8125419 0.8125419
## [1] 40.53391
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(NAcc_ISC) +
## scale(AIns_ISC) + scale(MPFC_ISC), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.38292 -0.39035 0.00291 0.42719 1.30755
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.4615011 0.4879555 25.538 3.83e-11 ***
## Theaters_US_W1_num 0.0011617 0.0001602 7.250 1.64e-05 ***
## scale(NAcc_ISC) 0.1735213 0.2448084 0.709 0.493
## scale(AIns_ISC) -0.0210687 0.2575483 -0.082 0.936
## scale(MPFC_ISC) -0.2471274 0.2289210 -1.080 0.303
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7769 on 11 degrees of freedom
## Multiple R-squared: 0.8362, Adjusted R-squared: 0.7766
## F-statistic: 14.04 on 4 and 11 DF, p-value: 0.0002672
## R2m R2c
## [1,] 0.7891992 0.7891992
## [1] 43.33375
# FW_US_M3: ISC data + affective data + behavioral data
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_ISC) + scale(AIns_ISC) +
## scale(MPFC_ISC), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.9034 -0.2448 -0.1094 0.2894 0.8893
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.2654495 0.3976897 30.842 1.94e-10 ***
## Theaters_US_W1_num 0.0012319 0.0001317 9.353 6.23e-06 ***
## scale(Pos_arousal_scaled) -0.2113917 0.2221020 -0.952 0.3661
## scale(Neg_arousal_scaled) 0.5028299 0.1826932 2.752 0.0224 *
## scale(NAcc_ISC) 0.4828712 0.2569446 1.879 0.0929 .
## scale(AIns_ISC) -0.0425640 0.2477111 -0.172 0.8674
## scale(MPFC_ISC) -0.4886649 0.2047304 -2.387 0.0408 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6035 on 9 degrees of freedom
## Multiple R-squared: 0.9191, Adjusted R-squared: 0.8652
## F-statistic: 17.05 on 6 and 9 DF, p-value: 0.0001896
## R2m R2c
## [1,] 0.8721157 0.8721157
## [1] 36.04087
# FW_US_M4: Whole data alone
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(NAcc_whole) +
## scale(AIns_whole) + scale(MPFC_whole), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.2756 -0.4046 0.2554 0.3573 0.8457
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.263e+01 4.182e-01 30.189 6.23e-12 ***
## Theaters_US_W1_num 1.103e-03 1.366e-04 8.072 6.00e-06 ***
## scale(NAcc_whole) 1.049e-01 1.960e-01 0.535 0.6030
## scale(AIns_whole) 6.497e-02 1.907e-01 0.341 0.7397
## scale(MPFC_whole) 3.494e-01 1.834e-01 1.906 0.0832 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6833 on 11 degrees of freedom
## Multiple R-squared: 0.8733, Adjusted R-squared: 0.8272
## F-statistic: 18.96 on 4 and 11 DF, p-value: 6.743e-05
## R2m R2c
## [1,] 0.8348434 0.8348434
## [1] 39.22441
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_whole) + scale(AIns_whole) +
## scale(MPFC_whole), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.72409 -0.20077 0.05612 0.32143 0.56474
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.5236381 0.3368316 37.181 3.65e-11 ***
## Theaters_US_W1_num 0.0011395 0.0001108 10.285 2.83e-06 ***
## scale(Pos_arousal_scaled) -0.0240594 0.1541740 -0.156 0.8794
## scale(Neg_arousal_scaled) 0.5143768 0.1705727 3.016 0.0146 *
## scale(NAcc_whole) 0.3995877 0.1855911 2.153 0.0597 .
## scale(AIns_whole) -0.1882468 0.1704791 -1.104 0.2981
## scale(MPFC_whole) 0.3402472 0.1519678 2.239 0.0519 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5317 on 9 degrees of freedom
## Multiple R-squared: 0.9372, Adjusted R-squared: 0.8954
## F-statistic: 22.4 on 6 and 9 DF, p-value: 6.255e-05
## R2m R2c
## [1,] 0.8995921 0.8995921
## [1] 31.98644
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(NAcc_onset) +
## scale(AIns_onset) + scale(MPFC_onset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.58549 -0.31399 0.00325 0.30015 1.33657
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.5939630 0.5219388 24.129 7.08e-11 ***
## Theaters_US_W1_num 0.0011143 0.0001729 6.443 4.79e-05 ***
## scale(NAcc_onset) -0.0383730 0.2490668 -0.154 0.880
## scale(AIns_onset) -0.2714296 0.3475997 -0.781 0.451
## scale(MPFC_onset) 0.1635024 0.3274233 0.499 0.627
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7902 on 11 degrees of freedom
## Multiple R-squared: 0.8306, Adjusted R-squared: 0.7689
## F-statistic: 13.48 on 4 and 11 DF, p-value: 0.0003201
## R2m R2c
## [1,] 0.782359 0.782359
## [1] 43.87564
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_onset) + scale(AIns_onset) +
## scale(MPFC_onset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.00596 -0.42891 0.06107 0.33857 1.20079
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.5037757 0.5296610 23.607 2.09e-09 ***
## Theaters_US_W1_num 0.0011466 0.0001757 6.525 0.000108 ***
## scale(Pos_arousal_scaled) -0.1620195 0.2546050 -0.636 0.540388
## scale(Neg_arousal_scaled) 0.3003505 0.2375941 1.264 0.237939
## scale(NAcc_onset) 0.2130005 0.3290135 0.647 0.533542
## scale(AIns_onset) -0.3247931 0.3756072 -0.865 0.409651
## scale(MPFC_onset) 0.0908062 0.3453251 0.263 0.798503
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7955 on 9 degrees of freedom
## Multiple R-squared: 0.8595, Adjusted R-squared: 0.7658
## F-statistic: 9.176 on 6 and 9 DF, p-value: 0.002047
## R2m R2c
## [1,] 0.7858827 0.7858827
## [1] 44.87972
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(NAcc_middle) +
## scale(AIns_middle) + scale(MPFC_middle), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.28065 -0.08859 0.02742 0.45957 0.71604
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.5103123 0.4667586 26.803 2.27e-11 ***
## Theaters_US_W1_num 0.0011442 0.0001542 7.419 1.33e-05 ***
## scale(NAcc_middle) 0.3449983 0.2323279 1.485 0.166
## scale(AIns_middle) 0.1234687 0.2221701 0.556 0.590
## scale(MPFC_middle) -0.1757236 0.2192511 -0.801 0.440
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7179 on 11 degrees of freedom
## Multiple R-squared: 0.8601, Adjusted R-squared: 0.8093
## F-statistic: 16.91 on 4 and 11 DF, p-value: 0.0001147
## R2m R2c
## [1,] 0.8184974 0.8184974
## [1] 40.80755
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_middle) + scale(AIns_middle) +
## scale(MPFC_middle), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.0602 -0.2823 0.1308 0.3278 0.7616
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.4515030 0.4946438 25.173 1.18e-09 ***
## Theaters_US_W1_num 0.0011653 0.0001639 7.109 5.61e-05 ***
## scale(Pos_arousal_scaled) -0.1373489 0.2486434 -0.552 0.594
## scale(Neg_arousal_scaled) 0.1460506 0.2456896 0.594 0.567
## scale(NAcc_middle) 0.2993163 0.2459860 1.217 0.255
## scale(AIns_middle) 0.1336529 0.3170095 0.422 0.683
## scale(MPFC_middle) -0.1413418 0.2571878 -0.550 0.596
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.7481 on 9 degrees of freedom
## Multiple R-squared: 0.8758, Adjusted R-squared: 0.7929
## F-statistic: 10.57 on 6 and 9 DF, p-value: 0.001212
## R2m R2c
## [1,] 0.8087687 0.8087687
## [1] 42.91168
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(NAcc_offset) +
## scale(AIns_offset) + scale(MPFC_offset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.81940 -0.37776 0.01451 0.25168 0.86566
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.3517604 0.3624849 34.075 1.67e-12 ***
## Theaters_US_W1_num 0.0012010 0.0001193 10.070 6.89e-07 ***
## scale(NAcc_offset) -0.1513350 0.1916306 -0.790 0.4464
## scale(AIns_offset) 0.1405807 0.1655208 0.849 0.4138
## scale(MPFC_offset) 0.5445114 0.1878277 2.899 0.0145 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5712 on 11 degrees of freedom
## Multiple R-squared: 0.9115, Adjusted R-squared: 0.8793
## F-statistic: 28.31 on 4 and 11 DF, p-value: 9.732e-06
## R2m R2c
## [1,] 0.8830375 0.8830375
## [1] 33.4901
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_offset) + scale(AIns_offset) +
## scale(MPFC_offset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.77884 -0.18659 0.05761 0.19377 0.57414
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.241e+01 2.958e-01 41.941 1.24e-11 ***
## Theaters_US_W1_num 1.182e-03 9.747e-05 12.126 7.05e-07 ***
## scale(Pos_arousal_scaled) 1.547e-01 1.331e-01 1.162 0.27513
## scale(Neg_arousal_scaled) 3.532e-01 1.307e-01 2.702 0.02433 *
## scale(NAcc_offset) -5.731e-02 1.588e-01 -0.361 0.72646
## scale(AIns_offset) 1.586e-02 1.425e-01 0.111 0.91384
## scale(MPFC_offset) 6.169e-01 1.570e-01 3.928 0.00347 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4622 on 9 degrees of freedom
## Multiple R-squared: 0.9526, Adjusted R-squared: 0.921
## F-statistic: 30.13 on 6 and 9 DF, p-value: 1.819e-05
## R2m R2c
## [1,] 0.9233876 0.9233876
## [1] 27.50086
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(NAcc_onset) +
## scale(AIns_middle) + scale(MPFC_offset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.96161 -0.26052 0.04203 0.26418 0.85995
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.256e+01 3.457e-01 36.343 8.24e-13 ***
## Theaters_US_W1_num 1.125e-03 1.138e-04 9.890 8.26e-07 ***
## scale(NAcc_onset) 7.063e-02 1.534e-01 0.460 0.65417
## scale(AIns_middle) 2.174e-01 1.467e-01 1.482 0.16644
## scale(MPFC_offset) 5.133e-01 1.526e-01 3.363 0.00633 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.544 on 11 degrees of freedom
## Multiple R-squared: 0.9197, Adjusted R-squared: 0.8905
## F-statistic: 31.49 on 4 and 11 DF, p-value: 5.737e-06
## R2m R2c
## [1,] 0.8935894 0.8935894
## [1] 31.93091
##
## Call:
## lm(formula = log(Gross_US_W1_num) ~ +Theaters_US_W1_num + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_onset) + scale(AIns_middle) +
## scale(MPFC_offset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.55141 -0.11755 -0.01153 0.16097 0.56232
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.243e+01 2.378e-01 52.247 1.73e-12 ***
## Theaters_US_W1_num 1.174e-03 7.844e-05 14.973 1.15e-07 ***
## scale(Pos_arousal_scaled) 7.111e-02 1.327e-01 0.536 0.605076
## scale(Neg_arousal_scaled) 4.784e-01 1.250e-01 3.827 0.004043 **
## scale(NAcc_onset) 2.905e-01 1.268e-01 2.292 0.047642 *
## scale(AIns_middle) 9.200e-03 1.297e-01 0.071 0.945004
## scale(MPFC_offset) 6.808e-01 1.170e-01 5.817 0.000254 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3699 on 9 degrees of freedom
## Multiple R-squared: 0.9696, Adjusted R-squared: 0.9494
## F-statistic: 47.89 on 6 and 9 DF, p-value: 2.521e-06
## R2m R2c
## [1,] 0.9503822 0.9503822
## [1] 20.37305
# Neuroforecasting: First Month US.
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.4772 -0.2704 -0.1041 0.2875 0.6834
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 15.1927135 0.2659405 57.128 5.46e-16 ***
## Theaters_US_M1 0.0001917 0.0000259 7.401 8.26e-06 ***
## scale(Pos_arousal_scaled) 0.1061626 0.1077262 0.985 0.344
## scale(Neg_arousal_scaled) -0.1161456 0.1040014 -1.117 0.286
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.3954 on 12 degrees of freedom
## Multiple R-squared: 0.8588, Adjusted R-squared: 0.8235
## F-statistic: 24.33 on 3 and 12 DF, p-value: 2.18e-05
## R2m R2c
## [1,] 0.8294972 0.8294972
## [1] 21.10915
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(NAcc_ISC) +
## scale(AIns_ISC) + scale(MPFC_ISC), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.4831 -0.2525 -0.1129 0.2019 0.8080
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.499e+01 2.894e-01 51.799 1.71e-14 ***
## Theaters_US_M1 2.131e-04 2.824e-05 7.548 1.13e-05 ***
## scale(NAcc_ISC) 1.574e-01 1.362e-01 1.155 0.272
## scale(AIns_ISC) -1.157e-01 1.406e-01 -0.823 0.428
## scale(MPFC_ISC) -3.612e-02 1.292e-01 -0.280 0.785
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.425 on 11 degrees of freedom
## Multiple R-squared: 0.8504, Adjusted R-squared: 0.796
## F-statistic: 15.63 on 4 and 11 DF, p-value: 0.0001645
## R2m R2c
## [1,] 0.8065301 0.8065301
## [1] 24.03129
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_ISC) + scale(AIns_ISC) +
## scale(MPFC_ISC), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.4909 -0.2308 -0.1194 0.3020 0.6478
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.513e+01 3.792e-01 39.891 1.95e-11 ***
## Theaters_US_M1 1.985e-04 3.802e-05 5.221 0.000549 ***
## scale(Pos_arousal_scaled) 4.136e-02 1.877e-01 0.220 0.830562
## scale(Neg_arousal_scaled) -1.348e-01 1.430e-01 -0.942 0.370562
## scale(NAcc_ISC) 7.412e-02 2.233e-01 0.332 0.747526
## scale(AIns_ISC) -1.100e-01 1.929e-01 -0.570 0.582576
## scale(MPFC_ISC) 4.567e-02 1.614e-01 0.283 0.783648
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4477 on 9 degrees of freedom
## Multiple R-squared: 0.8642, Adjusted R-squared: 0.7736
## F-statistic: 9.544 on 6 and 9 DF, p-value: 0.001772
## R2m R2c
## [1,] 0.7924233 0.7924233
## [1] 26.486
# FM_US_M17: Whole data alone
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(NAcc_whole) +
## scale(AIns_whole) + scale(MPFC_whole), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.48416 -0.25762 -0.05023 0.17241 0.66448
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.508e+01 2.883e-01 52.323 1.53e-14 ***
## Theaters_US_M1 2.031e-04 2.808e-05 7.230 1.69e-05 ***
## scale(NAcc_whole) 8.483e-02 1.265e-01 0.671 0.516
## scale(AIns_whole) 6.667e-02 1.261e-01 0.529 0.607
## scale(MPFC_whole) -1.177e-02 1.142e-01 -0.103 0.920
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4282 on 11 degrees of freedom
## Multiple R-squared: 0.8481, Adjusted R-squared: 0.7929
## F-statistic: 15.36 on 4 and 11 DF, p-value: 0.0001783
## R2m R2c
## [1,] 0.8037433 0.8037433
## [1] 24.27266
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_whole) + scale(AIns_whole) +
## scale(MPFC_whole), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.4385 -0.2695 -0.0189 0.2233 0.5891
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.544e+01 3.333e-01 46.326 5.09e-12 ***
## Theaters_US_M1 1.656e-04 3.334e-05 4.966 0.000774 ***
## scale(Pos_arousal_scaled) 1.405e-01 1.251e-01 1.123 0.290457
## scale(Neg_arousal_scaled) -2.440e-01 1.479e-01 -1.649 0.133464
## scale(NAcc_whole) -1.344e-01 1.674e-01 -0.803 0.442646
## scale(AIns_whole) 2.397e-01 1.535e-01 1.561 0.152892
## scale(MPFC_whole) 4.299e-02 1.158e-01 0.371 0.719043
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4024 on 9 degrees of freedom
## Multiple R-squared: 0.8903, Adjusted R-squared: 0.8172
## F-statistic: 12.17 on 6 and 9 DF, p-value: 0.0007106
## R2m R2c
## [1,] 0.8296253 0.8296253
## [1] 23.06842
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(NAcc_onset) +
## scale(AIns_onset) + scale(MPFC_onset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.43418 -0.31495 -0.04144 0.25868 0.66861
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.509e+01 3.108e-01 48.566 3.46e-14 ***
## Theaters_US_M1 2.020e-04 3.077e-05 6.564 4.06e-05 ***
## scale(NAcc_onset) 1.604e-01 1.328e-01 1.208 0.252
## scale(AIns_onset) 2.949e-02 1.747e-01 0.169 0.869
## scale(MPFC_onset) -1.160e-01 1.717e-01 -0.676 0.513
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4117 on 11 degrees of freedom
## Multiple R-squared: 0.8596, Adjusted R-squared: 0.8086
## F-statistic: 16.84 on 4 and 11 DF, p-value: 0.0001171
## R2m R2c
## [1,] 0.8178498 0.8178498
## [1] 23.01581
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_onset) + scale(AIns_onset) +
## scale(MPFC_onset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.46535 -0.28851 -0.09242 0.27231 0.55967
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.516e+01 3.447e-01 43.975 8.12e-12 ***
## Theaters_US_M1 1.954e-04 3.425e-05 5.703 0.000293 ***
## scale(Pos_arousal_scaled) 8.363e-02 1.440e-01 0.581 0.575569
## scale(Neg_arousal_scaled) -6.981e-02 1.338e-01 -0.522 0.614409
## scale(NAcc_onset) 7.689e-02 1.829e-01 0.420 0.683984
## scale(AIns_onset) 6.685e-02 2.022e-01 0.331 0.748547
## scale(MPFC_onset) -1.022e-01 1.934e-01 -0.529 0.609896
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4422 on 9 degrees of freedom
## Multiple R-squared: 0.8675, Adjusted R-squared: 0.7792
## F-statistic: 9.823 on 6 and 9 DF, p-value: 0.001594
## R2m R2c
## [1,] 0.7971308 0.7971308
## [1] 26.08615
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(NAcc_middle) +
## scale(AIns_middle) + scale(MPFC_middle), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.55081 -0.27048 -0.06784 0.26948 0.60581
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.529e+01 3.277e-01 46.660 5.36e-14 ***
## Theaters_US_M1 1.813e-04 3.257e-05 5.568 0.000168 ***
## scale(NAcc_middle) -1.685e-01 1.570e-01 -1.073 0.306087
## scale(AIns_middle) 9.364e-02 1.372e-01 0.683 0.509000
## scale(MPFC_middle) 1.430e-01 1.338e-01 1.069 0.307996
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4204 on 11 degrees of freedom
## Multiple R-squared: 0.8537, Adjusted R-squared: 0.8005
## F-statistic: 16.04 on 4 and 11 DF, p-value: 0.0001462
## R2m R2c
## [1,] 0.8105347 0.8105347
## [1] 23.67864
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_middle) + scale(AIns_middle) +
## scale(MPFC_middle), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.44224 -0.25199 -0.03754 0.31736 0.56953
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.541e+01 3.394e-01 45.411 6.09e-12 ***
## Theaters_US_M1 1.685e-04 3.387e-05 4.976 0.000764 ***
## scale(Pos_arousal_scaled) 7.843e-02 1.393e-01 0.563 0.587220
## scale(Neg_arousal_scaled) -1.364e-01 1.393e-01 -0.979 0.353147
## scale(NAcc_middle) -1.611e-01 1.572e-01 -1.025 0.332322
## scale(AIns_middle) 1.279e-01 1.829e-01 0.699 0.502230
## scale(MPFC_middle) 1.172e-01 1.466e-01 0.799 0.444654
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4199 on 9 degrees of freedom
## Multiple R-squared: 0.8805, Adjusted R-squared: 0.8008
## F-statistic: 11.05 on 6 and 9 DF, p-value: 0.001026
## R2m R2c
## [1,] 0.8155404 0.8155404
## [1] 24.43633
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(NAcc_offset) +
## scale(AIns_offset) + scale(MPFC_offset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.45430 -0.33861 -0.01704 0.20012 0.67277
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 15.0309917 0.3345914 44.923 8.13e-14 ***
## Theaters_US_M1 0.0002087 0.0000331 6.305 5.80e-05 ***
## scale(NAcc_offset) 0.0093962 0.1528761 0.061 0.952
## scale(AIns_offset) 0.0687571 0.1272437 0.540 0.600
## scale(MPFC_offset) -0.0271473 0.1510495 -0.180 0.861
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4457 on 11 degrees of freedom
## Multiple R-squared: 0.8355, Adjusted R-squared: 0.7757
## F-statistic: 13.96 on 4 and 11 DF, p-value: 0.0002736
## R2m R2c
## [1,] 0.7883149 0.7883149
## [1] 25.5531
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_offset) + scale(AIns_offset) +
## scale(MPFC_offset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.5148 -0.2385 -0.1090 0.2637 0.5562
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.516e+01 3.359e-01 45.117 6.46e-12 ***
## Theaters_US_M1 1.956e-04 3.335e-05 5.865 0.000239 ***
## scale(Pos_arousal_scaled) 1.061e-01 1.245e-01 0.852 0.416113
## scale(Neg_arousal_scaled) -1.534e-01 1.259e-01 -1.218 0.254159
## scale(NAcc_offset) -2.488e-02 1.539e-01 -0.162 0.875134
## scale(AIns_offset) 1.245e-01 1.303e-01 0.956 0.364142
## scale(MPFC_offset) -2.609e-02 1.520e-01 -0.172 0.867506
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4348 on 9 degrees of freedom
## Multiple R-squared: 0.8719, Adjusted R-squared: 0.7865
## F-statistic: 10.21 on 6 and 9 DF, p-value: 0.00138
## R2m R2c
## [1,] 0.8033241 0.8033241
## [1] 25.547
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(NAcc_onset) +
## scale(AIns_middle) + scale(MPFC_offset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.44882 -0.31788 -0.04391 0.21377 0.62567
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.524e+01 3.212e-01 47.436 4.48e-14 ***
## Theaters_US_M1 1.871e-04 3.189e-05 5.866 0.000108 ***
## scale(NAcc_onset) 1.630e-01 1.203e-01 1.355 0.202537
## scale(AIns_middle) 6.878e-02 1.147e-01 0.600 0.560913
## scale(MPFC_offset) 2.076e-02 1.265e-01 0.164 0.872627
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4149 on 11 degrees of freedom
## Multiple R-squared: 0.8574, Adjusted R-squared: 0.8056
## F-statistic: 16.54 on 4 and 11 DF, p-value: 0.0001272
## R2m R2c
## [1,] 0.8151494 0.8151494
## [1] 23.26338
##
## Call:
## lm(formula = log(Gross_US_M1) ~ +Theaters_US_M1 + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_onset) + scale(AIns_middle) +
## scale(MPFC_offset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.58968 -0.21370 -0.00677 0.24220 0.49517
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.536e+01 3.646e-01 42.125 1.19e-11 ***
## Theaters_US_M1 1.743e-04 3.652e-05 4.774 0.00101 **
## scale(Pos_arousal_scaled) 5.837e-03 1.593e-01 0.037 0.97156
## scale(Neg_arousal_scaled) -1.538e-01 1.573e-01 -0.978 0.35376
## scale(NAcc_onset) 9.125e-02 1.494e-01 0.611 0.55639
## scale(AIns_middle) 1.335e-01 1.670e-01 0.800 0.44441
## scale(MPFC_offset) -4.564e-02 1.600e-01 -0.285 0.78194
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4336 on 9 degrees of freedom
## Multiple R-squared: 0.8726, Adjusted R-squared: 0.7876
## F-statistic: 10.27 on 6 and 9 DF, p-value: 0.00135
## R2m R2c
## [1,] 0.8042702 0.8042702
## [1] 25.46328
# Total_US_M27: Seq data + affective data + behavioral data
##
## Call:
## lm(formula = log(Gross_Total_US) ~ +Total_Theaters_US + scale(Pos_arousal_scaled) +
## scale(Neg_arousal_scaled) + scale(NAcc_onset) + scale(AIns_middle) +
## scale(MPFC_offset), data = AllSubs_NeuralActivation)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.37432 -0.21825 -0.00611 0.20211 0.47060
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.568e+01 2.315e-01 67.724 1.69e-13 ***
## Total_Theaters_US 9.980e-05 1.411e-05 7.073 5.84e-05 ***
## scale(Pos_arousal_scaled) -3.391e-01 1.313e-01 -2.583 0.0296 *
## scale(Neg_arousal_scaled) 9.721e-02 1.518e-01 0.640 0.5380
## scale(NAcc_onset) 3.500e-01 1.273e-01 2.749 0.0225 *
## scale(AIns_middle) 1.651e-01 1.268e-01 1.302 0.2252
## scale(MPFC_offset) 2.330e-01 1.569e-01 1.485 0.1717
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.351 on 9 degrees of freedom
## Multiple R-squared: 0.9336, Adjusted R-squared: 0.8893
## F-statistic: 21.09 on 6 and 9 DF, p-value: 8.014e-05
## R2m R2c
## [1,] 0.8940089 0.8940089
## [1] 18.70016